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ABSTRACT 

We present Herschel PACS spectroscopy of the [01] 63 /im gas-line for three circumstellar disk sys- 
tems showing signs of significant disk evolution and/or planet formation: HR 8799, HD 377 and RX 
J1852.3-3700. [01] is undetected toward HR 8799 and HD 377 with 3a upper limits of 6.8 x 10" 18 
W m~ 2 and 9.9 x 10~ 18 W m~ 2 respectively. We find an [OI] detection for RX J1852. 3-3700 at 
12.3 ± 1.8 x 10~ 18 W vaT 2 . We use thermo-chemical disk models to model the gas emission, using 
constraints on the [OI] 63 /im and ancillary data to derive gas mass upper limits and constrain gas- 
to-dust ratios. For HD 377 and HR 8799, we find 3 a upper limits on the gas mass of 0.1-20 M®. For 
RX J1852. 3-3700, we find two distinct disk scenarios that could explain the detection of [OI] 63 /im 
and CO(2-l) upper limits reported from the literature: (i) a large disk with gas co-located with the 
dust (16-500 AU), resulting in a large tenuous disk with ~ 16 M® of gas, or (ii) an optically thick gas 
disk, truncated at ~ 70 AU, with a gas mass of 150 M®. We discuss the implications of these results 
for the formation and evolution of planets in these three systems. 

Subject headings: circumstellar matter — planetary systems: formation — stars: pre-main sequence 



1. INTRODUCTION 

Protoplanetary disks are found around most, if not all, 
forming s tars, with typical measured lifetimes of a few 
Myr (e.g. lHaisch et a"l.ll200l(: I Meyer et alj|2007j) . Planet 
formation is thought to be ongoing in these first few 
Myrs, leading to gas giant planets, ice-giants, and plan- 
etesimal belts that interact with the gas. The interaction 
with the gas feeds the growth of the planets, causes their 
migration, and circularizes their orbits. A key constraint 
governing the formation of gas giants, ice-giants and ter- 
restrial rocky planets, as well as their orbital evolution, is 
the timescale for gas clearing in a disk. Recently, the Ke- 
pler mission and radial velocity surveys have uncovered 
several candidate planets with sizes of order a few Earth 
that could be either r ocky "super Earths" or gaseous 
"min i Neptunes" (e.g. iFressin et al.l 120121 : iGaidos et al.l 
|20T1 [C7oll~et^[201ll) 

A sensitive search for the [SI] 25.23 /im and molecular 
hydrogen gas emission lines in the infrared with Spitzer 
of 10-100 Myr old sun-like stars with dusty disks has 
shown that less than ~ 13 - 190 M® of gas at 100 K was 
present, even for t he youngest stars in t he sample with 
ages of 5-20 Myr (Pascucc i et al.ll2006D . i.e. too little 
gas left to still form J upiter-like gas giants. Furthermore, 
iPascucci et al.l (j2006[ ) found CO derived gas limits of less 
than a few M® in the 10-40 AU region, for a sub-sample 
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of 8 disks with ages of 10 to a few 100 Myr, excluding 
even the further formation of ice-giants like Uranus or 
Neptune in these systems. More data on a larger sample 
are warranted to address the question of how frequent 
these young 10-100 Myr stars with debris disk systems 
still have enough remnant gas (e.g. ~10 M®) to form 
ice-giants like U ranus and Neptune- 
Work bv IGorti fe Hollenbachl (|2(M 120081k 
IHollenbach fc Gortil (I2009F has demonstrated that 
the [OI] 63 /im emission is one of the strongest gas 
emission lines originating from the 10-50 AU regions 
of protoplanetary disks, and a sensitive probe for 
measuring small amounts of remnant gas, tracing the 
gas disk as it dissipates. Numerical simulations show 
that major collisions probably o ccurred in the ic e -giant 
zone to 10 Myr o r longer (Grigorie va et al.l 120071 : 
iFord fc Chiang! 12007ft , beyond the lifetime of optically 
thick disks. Such activity would release additional 
gas in either a pseudo-continuous or episodic fashion. 
Observations of [OI] are needed to trace small amounts 
of gas, and help address questions such as: how long 
does gas capable of forming ice-giants and circularizing 
the orbits of forming terrestrial planets persist? Is gas 
generation in this zone strongly variable from one system 
to an other and hence possib l y epi s odic? ([Chen et all 
2001 lAlexander fc Armltagel \200% iBesla fc Wul 120071 : 
Kominami et a 1.1 120051) 



The PACS instrument (iPoglitsch et al.ll2010D on -board 
the Herschel Space Observatory (jPilbratt et all 12010ft 
provides a major opportunity to probe the far-IR wave- 
length regime that covers strong gas lines such as the [OI] 
63 /im line. The Herschel GASPS program is conducting 
a large scale survey of Herbig Ae, T Tauri, and debris 
disk systems, and has reported [OI] 63 /im gas detec- 
ti ons toward several disk sources with ages of ~ 4-12 Myr 
inlMathews et al.l (120101): iThi et al.l 1)201 If ) : IWoitke et al.l 
(|2011h : iTilling et all (|2012D . 

In this paper, we focus on three well-characterized disk 
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systems: HR 8799, HD 377, and RX J1852.3-3700. We 
describe our sample in Sect. 2, and the observations and 
data reductions in Sect. 3. In Sect. 4, we present thermo- 
chemical disk models and estimate gas mass (upper lim- 
its), and in Sect. 5, we discuss the implications for planet 
formation and evolution in these sources. 

2. SAMPLE 

For our sample selection, we focused on 10-100 Myr old 
debris disk systems that were well characterized with an- 
cillary observations and detailed dust disk models avail- 
able, and that would be well complementary to the ex- 
isting Herschel prog rams, in particular GASPS (e.g., 
iMathews et al.ll2010D . Each target either shows direct 
evidence of the presence of planets, or signs of significant 
evolution in the dust disk structure, such as asymmetries, 
gaps or holes, as derived from direct imaging and/or their 
SED. We summarize here our sample in stellar properties 
and signs of disk evolution / planet formation. 

HR 8799 — is a nearby (d=39.4 ± 1.1 pc. Ivan Leeuwenl 
2007) planetary system around an early- t ype star. We 
adopt the stellar properties derived by Idrav fc Kavd 
(1999). It is noted that while the effective tempera- 
ture is like that of an F0 star, it is observed to be a A 
Bootis type star, i.e. having anom alously weaker abun - 
dan ces in some of its metal lines, see lGrav fc Kavd (|1999T ) 
and lSadakand (|2006l ) for a description of the full spectral 
classification. We adopt the discussion of the age and 
potential associat ion with nearby stellar groups of HR 
8799 described in lHinz et al.l (|2010l ) and therefore adopt 
the age range estimate of 60^30° Myr from lMarois et ahl 
(2008). Four giant planets were discov ered through di- 
rect imaging ([Marois et al.l 120081 12010| ). with masses of 
9, 10, 10, and 7 Mj up , located at projec ted separations 
of 14. 5, 24, 38 and 68 AU respectively. iClose fc Males! 
(2010) searched for the presence of new wide companions 
and concluded there are no companions with masses >3 
M Jup within - 5" - 15" (~ 200 - 600 AU). 

HR8799 has a remnant debris disk, first ident i- 
fied with IRAS p hotometry (jZuckerman fc Sondl200l . 
I Chen et all (|2006| ) distinguished two components, a hot 
inner ring at 6 AU and a coo ler outer disk from Spitzer 
observations. ISu et al.l (|2009f ) revealed from Spitzer 70 
/jm images there to be a large halo of small dust grains 
at large radii, and they presented a disk model to con- 
sist of: 1) warm (T-150K) dust at 6-15 AU, within the 
innermost planet; 2) a broad 90-300 AU zone of colder 
dust (T ~ 45 K) beyond the outermost planet; and 3) a 
300-1000 AU halo of small grain s probably being blown 
out of the system (|Su et al.ll2009D . The disk was spatiall y 
resolved at 880 /xm observations bv lHughes et al.1 (|2011l ). 
and found to be consistent with the presence of the 2nd 
component of the Su model, albeit with a larger inner 
radius of 150 AU (rather than 90 AU). Pa tience et al.l 
(|2011h resolved the cooler debris disk at 350 fim, con- 
sistent with a 100-300 AU size, while also showing an 
indication of asymmetric structure. A dust mass for the 
debris disk of 0.1 M ffi has been derived based on sub- 
millimeter photometry ([Williams fc And rews 200j|). 

HD 377 — is a nearby (d=39.8 pc, Ivan Leeuwenl 12007) 
G2V star (~ 1.1 M Q ) with an age estimat e based on 
a variety of indicators, summarized in both lApai et al.l 



( 2008]) and lRoccatagliata et al.l (|2009l ). We adopt an age 
of 150t^ 5 (95% CL) Myr. 

HD 377 exhibits an SED consistent with a disk with 
a dust mass in small grains > 4 x 10~ 4 M 9 from < 10 
to > 100 AU (Hil lenbrand et al.l |20M ) and < 0.1 M m 
of wa rm gas (T > 100K) within 2 AU (|Pascucci et alJ 
120061). A detailed model including new sub-mm data 
(|Roccatagliata et al.l 120091 ) indicates 0.06 M e in grains 
< 3 mm located at orbital radii from 3-150 AU. It has 
not yet been resolved in scattered light. There is no 
evidence to date that it has formed a gas giant planet 
(|Apai et al.l 120081 ). 

RX J1852. 3-3700— is a K3 (1.1 M Q ) T Tauri star, in 
the gen eral vicinity of the CrA star-forming region (d = 
132 pc. iCarpenter et al.ll2005l ) - situated 1.84 deg away 
from the center of the CrA cluster, and co-moving with 
it in terms of both proper motio n and radial veloc i ty. W e 
adopt photometric colors from iNeuhauser et al.l {2000) , 
and estimate an average Av of 0.89 ± 0.05 mag. We 
adopt a spectral type of K3 with +- 1 subtype and derive 
a logT = 3.682 ± 0.020, and logi = -0.226 ± 0.087, 
wh ere L is in Lp, and T i s in K units. Comparison with 
the iBaraffe et al] (119981) . iD'Antona fc Mazzitellil (119971 ) 
and iSiess et ahl (jl997l ) evolutionary tracks gives us an 
range in age estimates of ~ 10 — 23 Myr. We adopt here 
an age of 18^ 2 Myr, derived fr om the Baraffc tracks. 

RX J1852. 3-3700 harbours a transitio nal disk with a 
low a ccretion rate (5 x 10 -10 M Q yr" 1 , iPascucci et al.l 
120071 ). In a comprehensive study utilizing the Spitzer 
SED, as well as sub-mm SMA and ASTE ob- 
servations, and detailed models based on RADMC 
(jDullemond fc Dominikl l200l . iHughes et ahl (|2010h in- 
ferred an optically-thin region within 16 AU and deter- 
mined that the gas content in the system is about an 
order of magnitude lower than expected for the standard 
assumed gas-to-dust mass ratio of 100 and a CO abun- 
dance of 10 -4 relative to H 2 , suggesting gas clearing is 
likely underway in this system. 

We summarize the key parameters of each source in 
Table [U 

3. HERSCHEL OBSERVATIONS 

Our targets were observed with Herschel PACS IFU 
spectrometer, with pointed chop/nod line spectroscopy 
AORs. The center spaxel (9.4" x 9.4") of the 5 x 5 array 
was centered on the source in both nod positions and we 
utilized a small chop throw (2"). The observation cen- 
tered on the [OI] 3 Pi 3 P 2 transition at 63.1852 ^m, 
exposing for 2 line repetitions and 5 cycles, resulting in 
a nominal on-source integration time of 3840 sec. Ob- 
serving date, program and AOR ID are listed in Table [2j 

We reduced all PACS spectra with the Herschel In - 
teractive Processing Environment (HIPE, lOttl l2010h . 
v7.0.0, user release, build 7.0.1931, with PACS calibra- 
tion tree version 26 Q. The spectra were extracted from 
the central spaxel, using the HIPE built-in correction for 
the point-source loss correction. 

At the time we carried out our spectral analysis, the 
latest HIPE version did not include formal propagation 

7 HIPE software obtained from 

http : //herschel . esac . esa. int/HIPE_download . shtml 
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Table 1 

Sample properties 



HR 8799 HD 377 RX J1852. 3-3700 



Spectral Type kA5 hFO mA5 v A Boo a G2V b K3 C 

Distance [pc] 39.4±l.l d 39.1 ± l.l d 130 c 

Age [Myr] ' 60t™ f 150l^ 5 E 18±3, 2 s 

Stellar Mass [M Q ] 1.47 a 1.02 h l.l c 

Stellar Temperature [K] 7430 a 5849 h 4808^ 

Stellar Luminosity [Lq] 4.92 a 1.17 h 0.59 g 

Log(L x / erg s" 1 ) [dcx] 28.11 1 29. l j 30.41 c 



Refe rences. — a . [Cray fc Kavc (1991), b. |M oorc & Paddock (1950), c. Neuhauscr ct al. (2000), d. van Lccuwcn (2007), c. 



Carpenter ct al. (2005), f. Marois ct al. (2008), g. this work, h. Casagrandc ct al. (2011), i. Robradc & Schmitt (2010), j. Wichmann ct al. 

(120031) 



Table 2 

Observation log 



Name 


Date Observed 


Program ID 


Observation ID 


HR 8799 


Jan. 1, 2011 


GTl_vgeers_l 


1342212242 


HD 377 


Jan. 11, 2011 


OTl_vgeers_2 


1342212530 


RX J1852. 3-3700 


Mar. 15, 2011 


OTl_vgeers_2 


1342216163 



of the flux errors on individual pixels. Therefore, we 
estimated the noise by measuring directly the dispersion 
of the pixels from a first-order polynomial that was fitted 
to the continuum over short spectral ranges around the 
expected central wavelength of 63.1852 /xrr0 of the [OI] 
emission feature (63.03-63.16 /im, and 63.22-63.39 /im). 
This pixel dispersion was checked for Gauss ianity with 
the Shapiro- Wilk test (jShapiro fc Willd 1 1965( 1 and found 
for all three sources to be consistent with being drawn 
from a normal distribution. The measured rms of the 
pixel dispersion is thus taken as a representative measure 
of the noise, and used further in our derivation of the line 
flux uncertainty and line flux upper limits. 

From the continuum fit, we report the continuum flux 
estimate or 3 a upper limits at 63.1852 /im in Table [3l 
We note the recommended PACS absolute calibration 
uncertainty for these observations is 30% but do not in- 
clude this in the uncertainties reported in Table [3l We 
note that for our one source with a significant detection, 
RX J1852. 3-3700, the estimated continuum flux at 63.19 
/im is consistent with the IRAS 60 /im and Spitzer MIPS 
70 /im measurements to within ~ 25%. 

To estimate the integrated line intensity of the detected 
[OI] feature, we first perform a least-squares fit of a 1st 
order polynomial to the continuum, and then remove this 
baseline continuum fit. Next, we use the same approach 
to fit a Gaussian to the remaining spectrum, restricting 
the central wavelength to 63.18-63.20 /im. Finally, the 
line flux is computed as the area of the best-fit Gaussian. 
The spectra, continuum and line fits are shown in Figure 
ffl 

The uncertainty on the line flux is evaluated using a 
Monte-Carlo approach. We start from the original spec- 
trum, which we assume to represent the mean, and then 
we add pixel-by-pixcl normally distributed random noise 
to the spectrum, where the amplitude of the random 
noise is set equal to the rms of the continuum flux dis- 
persion estimated previously. We then redo the line flux 

8 NIST Atomic Spectra Database v4.1.0 IIRalchenko et al.|[20TlT) 
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HR 8799 



cont. 1(7 rms: 0.133 Jy 



[01] 



- 



HD 377 



cont. lo rms: 0.122 Jy' 




RXJ 1852.3-3700 



cont. 1(7 rms: 0.090 Jy 




63.0 63.1 63.2 63.3 63.4 

Wavelength (/xm) 

Figure 1. Herschel PACS spectra for HR 8799, HD 377, and 
RX J1852. 3-3700, with a linear continuum fit subtracted. The 
continuum fit is based on two short wavelength ranges on either 
side of the [OI] feature. The rest wavelength of the [OI] 3 Pi — > 3 P2 
transition at 63.1852 (im is indicated with a dashed line. The 
dotted lines indicate hypothetical [OI] feature with line flux equal 
to the 3 <t upper limit and width equal to that derived for the RX 
J1852. 3-3700 detection. 



estimate on this new "resampled" spectrum, following 
the same method as described above. We repeat this 
procedure 5000 times and so create a distribution of esti- 
mated line fluxes for a distribution of resampled spectra. 
We take the robustly estimated dispersion of this distri- 
bution of line fluxes (using IDL's "biweight_mean" rou- 
tine) to be our formal uncertainty in the line flux. This 
method takes into account both the flux uncertainty in 
our spectrum (estimated from the rms) as well as the un- 
certainties on the fits to the continuum and the line. An 
example of this line flux distribution is shown in Figure 
[2] for RX J1852. 3-3700. 
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Figure 2. Distribution of line fluxes derived 1000 times from "re- 
sampled" spectra for RX J1852. 3-3700. The sigma of this distri- 
bution, derived with biweight_mean, is used as the approximation 
of the uncertainty of the line flux. 



We consider the [01] emission feature to be detected 
when the measured line flux is larger than 3 times the 
estimated uncertainty in the line flux. For sources with a 
non-detection, we estimate an upper limit on the line flux 
using the following approach. We start with the original 
spectrum, and add to it a fake [OI] feature in the form 
of a Gaussian, centered on 63.185 /im, and with width 
equal to 0.0683 /im, and we choose a specific amplitude 
taken from a predefined range of values. The value for 
width is taken from the best-fit to the source RX J1852.3- 
3700 for which [01] is clearly detected. Based on the 
spectrum with fake [OI] feature thus created, we repeat 
the measurement of the line flux and its uncertainty as 
described above, and evaluate whether this fake feature 
would have been considered a detection. We repeat this 
latter for a range of increasing amplitudes for the fake 
feature, and search so for the minimum feature amplitude 
at which we would have "detected" the feature (3tr). We 
take the line flux corresponding to this minimum feature 
amplitude to be our upper limit on the [OI] line flux. 

The [OI] 63 /im line is not detected toward HR 8799 
and HD 377 in the central spaxel, nor in the other spax- 
els. The derived 3u upper limits are reported in Table 
1 

The [OI] line is detected toward RX J1852.3-3700 with 
a line flux of 12.3 x 10~ 17 W m -2 and a la uncertainty 
of 1.8 x 10 -17 W m -2 . We note here too this reported 
uncertainty does not include the absolute PACS calibra- 
tion uncertainty of 30%. The [OI] line appears spectrally 
unresolved, with a FWHM ~ 77 km s _1 . Its line center 
wavelength of 63.191 /jm appears slightly redshifted to 
the rest wavelength, however, we note that this small 
shift is within la of the absolute wavelength calibration 
uncertainty at A = 63 /im, due to nominal pointing un- 
certainty of order 2". Inspection of the off-center IFU 
spaxels show no significant indication of extended [OI] 
emission. 

9 based on Fig. 12 and Sect. 5.2 of the PACS Spectroscopy Per- 
formance and Calibration document of May 31, 2011 
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4. MODELING 

4.1. Model description 

We calculate gas line emission us ing detailed thermo 
chemical disk models as described in Gorti fc Hollenbac 
([200l 12001 : lGorti et all (poll and lrlollenbach fc Gort 
(2009). We summarize the main features of the models 
here for completeness. Disk structure models are com- 
puted from an input surface density distribution (S(r), 
where r is the disk radius) by assuming vertical hydro- 
static equilibrium. Vertical gas density and tempera- 
ture structure is obtained by thermal balance. The main 
heating mechanisms considered are heating due to colli- 
sions with warmer dust grains, grain photoelectric heat- 
ing by Far Ultraviolet (FUV; 6eV< hu < 13.6 eV) pho- 
tons, X-rays, exothermic chemical reactions and cosmic 
rays. Gas cools by collisions with cooler dust grains 
and line emission by various ions, atoms and molecules 
whose abundances are calculated simultaneously by solv- 
ing for disk chemistry. Our chemical network consists of 
85 species and ~ 600 reactions that include photoioniza- 
tion and photodissociation processes. Freezing on grains 
to form ices is not included explicitly, but its effect on 
line emission is estimated when the dust temperatures 
are low enough for freeze-out. Gas radiative transfer is 
calculated by using an escape probability formalism and 
a simple two-layer model is used for the dust radiative 
transfer. 



4.2. Input parameters 

The stellar spectrum (optical to X-ray wavelengths) 
is a primary input to the model and obtained from the 
literature for the disks in this study. Dust parameters 
such as the size distribution, composition, mass and ra- 
dial extent are constrained by fitting the dust spectral 
energy distribution (SED). Extensive prior modeling of 
dust SEDs has been carried out for all the three disks 
considered, and we avail ourselves of these results. As 
the vertical structure in our models is determined by the 
gas density and temperature, we verify the fit to the dust 
SED after arriving at a disk solution. The abundance 
of PAHs, which affec t gas temperature via grain photo- 
electric heating (e.g.. lBakes fc Ti clens 1994]), is unknown 
and we have set this to zero. For the disks around HR 
8799 and HD 377, the dust grain size distribution inferred 
from SED modeling (|Su et al.ll2009t IRocc atagliat a et al.l 
2009) shows a lack of sub-micron sized grains and this is 
consistent with an insignificant PAH population in these 
disks. However, the disk around RX J1852. 3-3700 may 
have a more primordial dust grain size distribution and 
PAHs may be present. We adopt interstellar clou d gas 
phase abundances from Savage & Scmbach (1996). The 
adopted stellar and dust properties are listed in Tables Q] 
and S] respectively. 

4.3. Model Results: Gas Content of Disks 

HR 8799 — We consider that gas may be spatially co- 
located with the dust, either in the inner disk (6 < r < 
15 AU) or in the outer disk (90 < r < 300 AU) and de- 
rive upper limits to the gas mass. Although there could 
be gas present in the 15 — 90 AU region, this is dynami- 
cally unlikely due to the presence of multiple planets at 
these radii. We separately constrain gas mass in the in- 
ner and outer disks, i.e., we assume that gas is present 
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Table 3 

Summary of [OI] observations: observed central wavelength, line flux, continuum flux, and line luminosity. The [OI] and continuum 
detections are reported with 1 a uncertainty, the non-detections are reported as 3 a upper limits, and the reported uncertainties do not 

include the 30% absolute PACS calibration uncertainty, see Sect. [3] 



Name 


A 


Line Flux 


Cont. Flux 


Line Luminosity 




H 


[10- 18 W m- 2 ] 


[Jy] 


[1O" 6 L ] 


HR 8799 




< 6.8 


< 0.40 


< 0.32 


HD 377 




< 9.9 


< 0.37 


< 0.47 


RX J1852.3-3700 


63.191 


12.3 ± 1.8 


1.82 ± 0.09 


6.7 ± 0.4 



Table 4 

Adopted Dust Parameters 



Source 


Surface Density 


Radial Extent 


Dust size 


Dust Mass 


Reference 




E(r) <x 


[AU] 


H 


[M e ] 




HR 8799 












Inner Disk 


r° 


6 < r < 15 


1.5 < a < 4.5 


1.1E-6 


Su et al. f2009~) 


Outer Disk 


r° 


90 < r < 300 


10 < a < 1000 


0.12 


Su et al. f2009~) 


HD 377 


r° 


6.1 < r < 150 


14 < a < 3000 


0.058 


Roccatagliata et al. (2009) 


RX J1852.3-3700 


^180AU^ c -r/180AU 


16 < r < 500 


0.1 < a < 1000 


53.3 


Hughes et al. f2010) 



in either one of these regions. We compute the expected 
[OI] 63 /im line luminosity for disks by varying gas mass 
and thus derive upper limits based on the PACS non- 
detection for this source. 

We find that the gas mass has to be lower than ~ 
O.lMg in the inner disk and < 5 — 8M ffi in the outer 
disk in order to agree with the [OI] flux limits. For the 
inner disk model that corresponds to this upper limit, the 
densities and temperatures in the [OI] -emitting regions of 



the disk are 



10 8 cm~ J and T t 



gas 



100- 300 K. Heat 



released by photoreactions, mainly ionization of carbon, 
dominates gas heating. X-ray and FUV fluxes from the 
star are low and the absence of PAHs and low dust colli- 
sional cross-section per H atom lead to low efficiency in 
heating the gas. [OI] is the main coolant. In the outer 
disk model corresponding to the upper limit, the disk 
is extended and cold, with the gas temperature in the 
[OI] emitting region T gas ~ 40 — 70 K, and gas densities 
n ~ 10 6 cm -3 . Photoionization of carbon again domi- 
nates gas heating and [OI]63/im and [CII] 158 /im lines 
are strong coolants. The gas/dust ratio upper limit in 
the outer disk is consistent with the canonical interstel- 
lar value of 100, while we find no useful constraint on the 
inner disk (< 10 5 ). 

HD 377 — In addition to the limit on [OI] emission pre- 
sented in this paper, this source has flux limits on H2 
S(l) emission and CO(2 — 1) and CO(3 — 2) rotational 
line emission of 2.2 x 10~ 18 , 5 x lCT 20 a nd 2.2 x 10~ 19 W 
m~ 2 respectively ()Pascucci et al.l [20061) . We derive an 
upper limit on the gas mass between 6 — 150 AU of 
~ 1 M ffi based on all these constraints, and a correspond- 
ing gas/dust ratio limit of < 17. We note that the con- 
straint derived from [OI] is more stringent than the limit 
of < 20 M m for T=1 00 K gas based on H 2 flux limits 
(jPascucci et al.ll2006t) . 

RX 31852.3-3700— The PACS [OI]63Aim detection, 
Spitzer IRS [Nell] 12. 8 /xm detection, and the non- 
detection of both H2 S(l) and CO (2 — 1) emission are 
all used for comparison with model calculated line fluxes 



to determine the distribution of gas in the disk. The 
fluxes and upper limits for these lines are 1.2 x 10 _17 W 
m~ 2 for [OI] (this pap er), 4.3 x 10~ 18 W m" 2 for [Nell] 
(|Pascucci et all 120061), < 1.0 x 10~ 18 W m" 2 for H 2 
S(l) (IPascucci et al.l[200l and < 8 x 10~ 22 W m" 2 for 
COf2^lT (|Hughes et allf2010l) . 

First, we consider constraints set by [OI] and CO emis- 
sion and models where the inner hole region (r < 16 AU) 
is devoid of gas. [OI] and CO(2 — 1) emission arises from 
relatively cool gas and these lines are better at tracing 
the disk mass. There are then two possible disk config- 
urations that can give rise to the observed [OI] and are 
still consistent with the non-detection of CO. 

In one case, the gas disk is extended and spatially co- 
located with the dust (up to 500 AU) and the mass of gas 
is very low, with the [OI] emission being optically thin. 
In such low mass, tenuous disks, CO is photodissociated 
by stellar and interstellar UV photons and there is no 
detectable CO(2 — 1) emission. We derive a gas mass of 
~ 16 within 16 — 500 AU in this case, and a highly 
reduced gas/dust mass ratio < 1. We note that although 
this model disk extends to 500 AU, both the CO(2 - 1) 
and [OI] emission decline sharply beyond r ~ 200 AU and 
we do not probe gas in the 200 — 500 AU regions. Higher 
disk masses are possible, if we consider the possibility 
that CO might be frozen out on cold dust grains. In 
models that include the effects of photodesorption, and 
assume freezing of CO and of O as water ice when the 
dust temperature is below some critical value, we esti- 
mate slightly higher disk masses. We assume that CO 
and H2O ice freeze at dust grain temperatures of 20 K 
and 80 K respectively, and account for ph otodesorption 
of ices as described in iGorti et al.l (| 2 1 If ) . With these 
assumptions, we find a maximum disk mass of ~ 40 M®, 
and a gas/dust mass ratio of ~ 3 in the disk extending 
from 16-500 AU. 

An alternate model is a truncated gas disk with the 
fiducial gas/dust mass ratio of 100 and is optically thick. 
This model disk produces excess [OI] and CO(2 — 1) emis- 
sion if it is as large as the dust disk and has to be limited 
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to a radial extent of ~ 70 AU to explain the observa- 
tions. The gas mass for the best fit model is 150 
from 16 — 70 AU. Freezing effects are not significant for 
the warmer dust temperatures of the compact disk con- 
figuration. There is also a possibility that CO might be 
under-abu ndant in the disk, as is often inferred for other 
disks (e.g. jQi et aIll2006ll20TlT) . In such a case, we infer 
CO depletion factors of ~ 30 — 100 in the disk around 
RX J1852. 3-3700. We next include the constraints intro- 
duced by mid-infrared line emission. The non-detection 
of H 2 S(l) excludes the possibility of a "full" gas disk, 
i.e., a disk of mass 0.016 M (5.3 x 10 3 M®) with a 
gas/dust ratio of 100. All disk configurations with masses 
a factor of 10 or more lower are possible within the con- 
straints set by H2 S(l), corresponding to upper limits to 
th e gas/dust rat i o of < 10, consistent with the findings 
of IHuehes et~aTl (2010). [Nell] emission is not a sen- 
sitive tracer of disk mass, and is reproduced to within 
a factor of ~ 3 in all our disk models. Model calcula- 
tions always give lower [Nell] emission than observed in 
the models that assume a 16 AU hole, which suggests 
that trace amounts of gas may be present within the 
hole. iPascucci et al.l (|2007|) estimate an accretion rate of 
5 x 10 _10 M©yr~ 1 from the Ha equivalent width, indi- 
cating the presence of at least some gas at r < 16 AU. 

We calculate disk life times for the models 
using our time-dependent photoevaporatio n models 
(|Gorti fc Hollenbachl [2009t IGorti et~aLl [2001 ) . For the 
dense, truncated disk extending from 16 — 70 AU we find 
a disk lifetime of ^0.5 Myr, while for extended, tenuous 
low mass disk we find a very short disk lifetime of only 
~ 4 x 10 4 yr. This argues in favor of the dense truncated 
disk model as the more likely scenario. 

5. DISCUSSION 

Here we explore some of the implications of our find- 
ings. We first compare our [OI] detections and upper 
limits with what has so far been found for the GASPS 
sample. Then, we explore scenarios for potential further 
planet formation in our target sample. 

5.1. Comparison with the GASPS program 

The Herschel GASPS program is carrying out a gas 
line emission survey toward a sample of 240 young stars, 
focussing primarily on [OI] 63.185 /im, [OI] 145.525 /mi 
and [CII] 157.741' fxra. GASPS have presented [OI] 63 
/xm gas observ ations toward 7 disk s o urces with ages o f 
~ 4-12 Myr inlMathews et al.l (T2 010): T hi et al.l (|20H ; 
IWoitke et all (12011ft : iTilling et al.l (12012ft . [OI] 63 /im 
has been detected toward 6 of these 7 sources, with 1 
non-detection for HD 181327. For the [OI] 63 detections 
in GASPS, the reported lines fluxes are typically a factor 
of a few stronger than what we find for the ~ 18 Myr old 
T Tauri disk RX J1852.3-3700 (12.3 xl0~ 18 W m" 2 ). 
Of particular note is the 4 Myr old Herbig Ae star HD 
163296, which displays a very strong [OI] 63 /jm line 
flux of - 193 x 10" 18 W m" 2 . The [OI] 63 /im upper 
limits that we find for our two older (60-150 Myr) disk 
sources are comparable to the < 9.3 x 10~ 18 W m~ 2 
presented for HP 181 327 (12 Myr old F5 star at 50.6 pc) 
(|Mathews et alj|2010j) . 

[OI] 63 fim is detected in one of the transitional disks 
in the GASPS sample, TW Hya, a 10 Myr K6V star with 



an inner 4 AU cavity in its disk. iThi et al] (|2010ft find 
that while the disk is still massive (0.5-5 x 1O _3 M0), a 
significant fraction of the primordial gas has already dis- 
appeared (gas-to-dust ratio between 2 and 26), possibly 
due to e vaporation by the strong X-ray flux from TW 
Hya (cf. IPascucci fc Sterzikl [2009ft . 

For HD 169142 , a ~ 6 Myr Herbig Ae disk, 
IMeeus et all (|2010D also detect [OI] 63 /im gas, and they 
propose this disk may be transitional. With a possible 
5-20 AU hole, and a reduced gas-to-dust ratio between 
22 and 50, the disk is still quite gas rich with a gas mass 
of 3-7 x 10" 3 M Q . 

We note that the results from our [OI] 63 fj,m observa- 
tions are overall consistent with those reported from the 
GASPS survey so far. More data are needed to derive 
a significant empirical relation for the [OI] 63 /im line 
strength as a function of age and stellar mass. 

5.2. Implications for planet formation 

With strong limits on the remnant gas in our sample 
of 3 evolved circumstellar disks, we can explore the pos- 
sibilities for on-going planet formation. To address this, 
we estimate the size that we could expect proto-planets 
to grow through collisional accretion ("isolation mass") 
given the ages of these systems, at the radii that a de- 
bris disk is still observed to be present. To estimate the 
isolation mass as a funct ion of orbital rad ius, we use the 
standard prescription by llda fe; Lml (|2008ft : 

„ n1R 3/2,3/2/ a N3/4-(3/2)( 9d -3/2) / M \-V2 

(1) 

where we assume q c i = —1 for the slope of the power 
law for the dust mass surface density E<j, where fd is 
a multiplicative factor used to scale E^, and ?7i C o is a 
step function describing the change in dust surface mass 
density at the ice-line. From this, the isolation mass 
scales with orbital radius and stellar mass as 

M Ciiso oc a 4 - 5 M- 1/2 (2) 

Next, we estimate what is the timescale for growing a 
planet up to this isolation mass through growth of solids, 
for given orbital radius and stellar mass, with the relation 



from rGoldrcic h~et al] ([20040, with planet density p p , 
planet radius R p , surface mass density E^, and planet 

orbital frequency Cl p . With Q p oc a~ 3 / 2 All^ 2 , assuming 
E<j cx a~~ 1 M*, and setting the planet mass to the isolation 
mass, it follows that this timescale scales with orbital 
radius and stellar mass as: 

t cx a 4 M- 5/3 (4) 

We normalize these relations to the mass surface density 
of a disk that would be capable of forming an ~1 M® 
isolation mass core at 5 AU (beyond the ice line) in 1 
Myr around an 1 M star, and then scale these relations 
to the different stellar masses and disk radii appropriate 
for our sample. 

10 this approximation does not include gravitational focussing, 
which means these are rather upper limits on the timescale 
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Following these approximate estimates for isolation 
mass and formation timescales, and given the ages of 
our sample disks, we find the following: 

For HD 377, at an age of 150 Myr, several cores of a 
few earth mass could have formed already within the 6- 
150 AU debris disk. Our gas mass upper limit of < 1 M® 
suggests there is not enough gas mass left to still form 
any more ice giants, but our upper limits cannot rule 
out the possibility that some existing gas / ice giants or 
dwarf planets could still be accreting remnant gas. The 
H2 and He might accrete as a thin atmosphere, while the 
"metal" content (Z > 2) would be expected to freeze out 
as solids and/or form an icy surface layer. 

In HR 8799, with a derived gas mass upper limit of 
< 5 - 8 M® in the outer disk (90-300 AU), there is not 
enough gas still present beyond 90 AU for the further 
in-situ formation of gas giants. At the age of 60 Myr, 
only small dwarf protoplanets (of order 10 times smaller 
than Pluto size) could be expected to have grown through 
collisional accretion in the debris disk beyond a radius 
of 90 AU. With the upper limit of 5-8 M®, such small 
protoplanets would accrete some of this material to form 
a solid icy surface layer. 

In the RX J1852. 3-3700 system, we have distinguished 
two possible disk scenarios that could be consistent with 
the detected [OI] emission. In the dense truncated disk 
scenario, with an estimated gas mass of ~ 150 M® and 
gas/dust mass ratio of ~ 100, there could still be enough 
gas present to form icy giants or volatile rich super-earths 
beyond 16 AU. The isolation mass beyond 16 AU is es- 
timated to be of order a hundred Earth masses and this 
could be reached on timescales of ~ 100 Myr. At the es- 
timated age of ~18 Myr, a few M® core may have formed 
beyond 16 AU, for which there is probably enough gas 
left to yield a final gas to dust ratio of 1:1 as inferred 
for the ice giants in our solar system. In this dense disk 
scenario, the inner hole may have been carved out when 
the disk was still quite massive (throu gh planet formation 
and /or X-ray photo-eva poration, e.g. lMarsh fc Mahonevl 
Il992tl0wen et al.ll2011[ ). whil e the outer disk was ablated 
by FUV photo-evaporation (jGorti fc Hollenbachl l2009h . 
The gas would have been present in the giant and ice 
giant planet formation regions of a putative planetary 
system, even at very late stages, affecting gas content 
and dynamics. The gas disk would have dispersed by 
erosion at both the inner and outer edges and evolve 
into a torus-like structure like we see at present. In the 
alternative scenario, we considered an extended gas disk 
co-located with the dust out to 500 AU with a much lower 
gas mass of 16 M® and highly reduced gas/dust mass ra- 
tio (< 1). In such a disk, there would be too little gas 
left for the further formation of Jovian-type planets, al- 
though potentially enough to still influence the dynamics 
of already formed planets. The surface mass density of 
the disk could have dec lined through photo-ev aporation 
(e.g., by EUV photons. I Alexander et al.ir2006|) . and vis- 
cous spreading, growing the disk to large radii. 

For future investigation, it will be interesting to ex- 
plore the effect of gas drag on the particle growth as 
well as the damping of orbital eccentricities in these ten- 
uous outer disks, and to consider global scenarios for 
planetary architecture that take into account the self- 
consistent formation and evolution (i.e. migra t ion) s ce- 
narios as described in e.g. lWalsh fc Morbidellil (|2011| ). 



6. CONCLUSIONS 

We present Herschel observations of the [OI] 63 /im 
gas line for three older debris disks, HR 8799, HD 377, 
and RX J1852.3-3700. [OI] is undetected toward HR 
8799 and HD 377 with 3ct upper limits of 6.8 x 10" 18 
W m~ 2 and 9.9 x 10~ 18 W m~ 2 respectively, while 
detected toward RX J1852.3-3700 with a line flux of 
12.3 ± 1.8 x 10 -18 W mT 2 . We use thermo-chemical 
disk models to model the gas emission, using constraints 
on the [OI] 63 /im and available ancillary data on the 
dust structure, to derive gas mass upper limits and con- 
strain gas-to-dust ratios. Based on the non-detections 
for HR 8799 (60±£g° Myr) and HD 377 (150t£j 5 Myr), 
we derive gas mass upper limits between 0.1-1.5 M® for 
these sources. For the [OI] detection toward RX J1852.3- 
3700 (18 Myr), we infer two possible scenarios for either 
a smaller 16 — 70 AU optically thick disk of 150 M® or 
alternatively a large 16 — 500 AU optically thin disk with 
about ~ 16 M® gas mass that has a significantly reduced 
gas/dust ratio of ~1. It appears that there is no possi- 
bility for further ice giant formation for HD 377 and HR 
8799, which could only still be forming smaller dwarf 
planets. In RX J1852. 3-3700, there appears to be suffi- 
cient solid and gaseous material for an ice giant to still 
form beyond 16 AU. 
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